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In brief 

Wheatcroft et al. show that nestling pied 
flycatchers preferentially beg in response 
to songs from their own population’s 
dialect. Nestlings beg weakly to songs 
from other populations and from collared 
flycatchers, a closely related species. 
These findings suggest that cultural 
evolution drives the emergence of pre- 
mating reproductive barriers. 
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SUMMARY 


Oscine birds preferentially respond to certain sounds over others from an early age, which focuses subse- 
quent learning onto sexually relevant songs. Songs vary both across species and, due to cultural evolution, 
among populations of the same species. As a result, early song responses are expected to be shaped by se- 
lection both to avoid the fitness costs of cross-species learning’ and to promote learning of population- 
typical songs.? These sources of selection are not mutually exclusive but can result in distinct geographic 
patterns of song responses in juvenile birds: if the risks of interspecific mating are the main driver of early 
song discrimination, then discrimination should be strongest where closely related species co-occur.^ In 
contrast, if early discrimination primarily facilitates learning local songs, then it should be tuned to songs 
typical of the local dialect." Here, we experimentally assess the drivers of song discrimination in nestling 
pied flycatchers (Ficedula hypoleuca). We first demonstrate that early discrimination against the songs of 
the closely related collared flycatcher (F. albicollis) is not strongly affected by co-occurrence. Second, across 
six European populations, we show that nestlings' early song responses are tuned to their local song dialect 
and that responses to the songs of collared flycatchers are similarly weak as to those of other conspecific 
dialects. Taken together, these findings provide clear experimental support for the hypothesis that cultural 
evolution, in conjunction with associated learning predispositions, drives the emergence of pre-mating 
reproductive barriers. 


RESULTS AND DISCUSSION 


Species discrimination without species interactions 

In pied flycatchers, song is produced by males on their arrival to 
the breeding grounds, during territorial interactions and pair for- 
mation,? and plays an important role in female mate choice.?'? 
Pied flycatchers co-occur with closely related collared flycatchers 
(F. albicollis) in central Europe as well as the Baltic islands of 
Gotland and Oland.'' The two species are highly similar in their 
morphology, ecology, and plumage, and they occasionally 


hybridize.'' Hybrids are completely sterile, '? implying large costs 
to cross-species song learning.'^ Where they co-occur, pied 
flycatcher nestlings discriminate against the songs of collared fly- 
catchers, an ability that arises independently of the early environ- 
ment.'^ If early song discrimination is primarily influenced by se- 
lection to avoid cross-species mating interactions, nestlings 
should more strongly discriminate against collared flycatcher 
songs where the two species co-occur (Figure 1C). 

Since songs vary within and across populations of both spe- 
cies,'? we analyzed songs from 168 individual pied and collared 


Current Biology 32, 1-6, December 5, 2022 © 2022 The Author(s). Published by Elsevier Inc. 1 
This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). 


Please cite this article in press as: Wheatcroft et al., Species-specific song responses emerge as a by-product of tuning to the local dialect, Current 


Biology (2022), https://doi.org/10.1016/j.cub.2022.09.063 


© Cell?ress 


OPEN ACCESS 


0.0 


PC1 (19.5%) 


0.1 


Current Biology 


Figure 1. Study populations, song variation, 
and hypotheses 

(A) Study sites: England (UK), Spain-Valsaín (SV), 
Spain-La Hiruela (SL), the Netherlands (NL), Oland 
(OL), Tovetorp (TO), and Finland (FI). Yellow denotes 
breeding range of pied flycatchers, blue denotes 
breeding range of collared flycatchers, and green 
denotes where both species co-occur in sympatry. 
Starred green circles denote sympatric locations 
used to estimate each population's distance to 
sympatry. 

(B) Species and population variation in song. PC1 
and PC2 scores for song phrases, averaged within 
individual males (see Figure S1A). Colors refer to 
populations in (A). Gray denotes recordings ob- 
tained from additional populations. Large circles 
denote population means. Whiskers give SE around 
each population mean. Ellipses contain 95% of the 
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flycatcher males in order to obtain a measure of song variation 
independent from species category (Figure S14). Pied flycatcher 
song phrases consist of 2-2.5 s long sequences of around 9 
smaller units, termed syllables, that are often repeated within 
and across song phrases. ^'^ We took multiple acoustic mea- 
sures of the syllables in each individual's song and found that 
PC1, a single principal component summarizing acoustic varia- 
tion in song phrases, correctly classified 10096 of the collared 
flycatcher singers (N = 41 individuals, mean PC1 score + SD, 
0.11 + 0.02) and 97.6% of the pied flycatcher singers (N = 127, 
—0.02 + 0.04) to the species level (Figure 1B; https://doi.org/ 
10.17632/v49k6v4fw3.1). Given this independent support for 
classifying songs at the species level, we considered playback 
species as a categorical variable when explaining song re- 
sponses. We found strong evidence that pied flycatcher nes- 
tlings discriminate in favor of pied flycatcher songs. In an exten- 
sive playback experiment across six populations (Figure S1B; 
Video S1), nestlings produce more begging calls in response 
to playbacks of pied flycatcher songs (estimated coefficient on 
pied flycatcher playbacks: 0.33 + 0.10 SE, likelihood-ratio test: 
xX? = 10.73, p = 0.001; Figure 2; Table S1; https://doi.org/10. 
17632/v49k6v4fw3.1). We found little evidence for an effect of 
age on species discrimination (likelihood-ratio test: X? = 0.11, 
p = 0.74), and discrimination was similar across the tested nest- 
ling age categories (Table S1). 

Next, we evaluated the effect of co-occurrence on early spe- 
cies discrimination by comparing nestling responses in allop- 
atry with those previously collected in sympatry on Öland. We 
found no evidence of an effect of co-occurrence on species 
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individuals for each species. Large yellow and blue 
circles denote species means. 

(C) Overview of hypotheses explaining species 
discrimination in nestling songbirds. (I) The costs of 
interspecific mating in sympatry drive selection to 
reduce cross-species learning, which, in turn, leads 
to early species discrimination. (Il) Selection to learn 
population-typical songs promotes early dialect 
discrimination. Due to larger song differences be- 
tween species than among dialects, this leads to the 
pattern of species discrimination even in allopatric 
populations. Colors of boxes denote whether the 
hypothesis predicts species discrimination in sym- 
patry (green) or both allopatry and sympatry 
(mixed). 


-0.05 -0.10 


discrimination (likelihood-ratio test: X? = 0.48, p = 0.49). Post 
hoc tests suggest stronger responses to pied flycatcher play- 
backs in both allopatry and sympatry with collared flycatchers 
(Figure 2; Table S2). Although typical pied flycatcher natal 
dispersal distances are between 5-20 km,'®?? there are 
extreme reports of up to 660 km.?^' We thus considered that 
populations even outside of current zones of contact might 
be or have been in the recent past exposed to species interac- 
tions. Therefore, we assessed the impact of geographic dis- 
tance to sympatry on nestling song responses. We found little 
evidence for a positive effect of proximity to sympatry on spe- 
cies discrimination (likelihood-ratio test: X? = 0.60, p = 0.44): 
discrimination was similar at varying distances to sympatry 
(Table S2). 

These results run counter to the idea that early discrimination 
against collared flycatcher songs arises primarily due to 
selection against cross-species learning. Instead, our results 
demonstrate that the ability to discriminate collared from pied 
flycatcher songs arises independently of co-occurrence. 
Although species interactions are hypothesized to impact re- 
sponses to sexual traits in general,’ we have shown here 
that "species discrimination" in nestling birds is a pattern that 
has arisen in the apparent absence of them (i.e., *endoge- 
nously”).2° However, species interactions may potentially influ- 
ence discrimination both directly and indirectly in other ways. 
First, experience can modify the strength of responses to other 
species’ sexual signals.^^ On Öland, where both pied and 
collared flycatchers co-occur, adult female pied flycatchers ex- 
press stronger preferences for pied flycatcher songs compared 
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Figure 2. Song discrimination in allopatry and sympatry 

Nestling responses in allopatry (left) and sympatry (right) to playbacks of 
collared and pied flycatcher songs (see Figure S1B). Clutches that produced 
no begging calls during any playback were removed for plotting purposes. 
Results are averaged across nestlings within nest. Boxplots give median, 
quartiles, and quartiles + 1.5 x IQR. p values are from post hoc tests (see 
Tables S1 and S2). 


with in allopatry,'° and collared flycatcher males have more 
species-specific song responses compared with juveniles? 
(see also Hudson et al.^5), suggesting that experience gained 
throughout development may strengthen the song discrimina- 
tion we observe in nestling birds. Second, discrimination 
against the songs of collared flycatchers could have evolved 
in pied flycatcher nestlings as a by-product of selection against 
learning the songs of other species that co-occur with pied fly- 
catchers throughout their breeding range, such as the great tit. 
Even in this scenario, we have shown that interactions with 
collared flycatchers in particular are not required for the emer- 
gence of discrimination against their songs.^? 


Song discrimination tuned to the local dialect 
Adult birds are thought to gain social and sexual benefits by 
producing and preferring song types that are common in their 
population,9'^ which drive discrimination toward locally typical 
songs.^/ ^? However, social transmission of songs promotes cul- 
tural evolution within populations,7°"°° leading to the prediction 
that early song responses should themselves evolve to track 
ongoing song evolution.^/ Co-evolution between songs and 
associated learning predispositions might therefore promote 
discrimination against the songs of non-co-occurring species 
and even among populations of the same species.”° 

We evaluated the hypothesis that early song discrimination 
arises due to selection to learn locally common songs by testing 
the prediction that nestlings respond strongest to songs typical 
of their own population (Figure 1C). Song acoustic features 
vary statistically across pied flycatcher populations. The majority 
(82.696) of song phrases could be correctly assigned to popula- 
tion in a linear discriminant analysis (Figure 1B; https://doi.org/ 
10.17632/v49k6v4fw3.1), providing clear support for the exis- 
tence of song dialects. Next, we assessed whether nestling 
song responses are tuned to their own population's song dialect. 
We related nestling begging responses to a given playback 
based on the acoustic distance between the playback and the 
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nestlings’ local dialect ("local-similar," “local-dissimilar,” 
*collared-similar"; Figure 3A). We found a clear effect of similar- 
ity to the local dialect on responses (likelinood-ratio test: X = 
6.74, p = 0.034), with post hoc tests indicating that nestlings 
respond strongest to local-similar playpacks and weakest to 
local-dissimilar and collared-similar playbacks (Table S3). In 
addition, there was a clear effect of age on these responses (like- 
lihood-ratio test: X = 9.93, p = 0.007). Although there was no 
clear effect of age on discrimination between local-similar and 
collared-similar playbacks, discrimination between local-similar 
and local-dissimilar songs increased with nestling age 
(Table S3), with a clear difference between ages 11 and 12, the 
age categories with the largest sample sizes (interaction 
contrast: p = 0.004; Figures 3B and 3C). To account for the pos- 
sibility that increased responses to local-similar songs with age 
are a consequence of the large difference in sample sizes of nes- 
tlings aged 11 and 12, we bootstrapped our results using sub- 
sampling. Models run on 99 of 100 sub-sampled datasets 
showed significantly stronger discrimination between local- 
similar and local-dissimilar songs at age 12 compared with age 
11 (interaction contrasts: p « 0.015; https://doi.org/10.17632/ 
v49k6v4fw3.1), implying that differences in song responses 
across age are not due to variation in sample size. 

A long-standing debate has centered around the role of song 
dialects in the emergence of incipient reproductive barriers.?' 
Some recent studies have suggested that dialects may influence 
patterns of mating in mixed populations. In white-crowned spar- 
rows, dialect divergence is in part linked with genetic diver- 
gence, and song acts as a potential behavioral barrier to repro- 
duction between subspecies.?^ Moreover, juvenile sparrows 
preferentially learn their own dialect,?? suggesting direct links 
between divergence in songs and early discrimination. We eval- 
uated the association between dialect divergence, genetic diver- 
gence, and geographic distance using previously published pair- 
wise Fs measures for five populations included in our study.?* 
We found support for a model including both genetic and 
geographic distance (r° = 0.60, F = 5.18, p = 0.01), which demon- 
strated a clear positive association between dialect and genetic 
divergence (t = 3.11, p = 0.010) and a trend for a negative asso- 
ciation with geographic distance (t = —2.17, p = 0.088). Since we 
show that early nestling responses to songs are tuned to their 
local dialect, these results are consistent with the idea that ge- 
netic variation among populations might directly impact both 
song dialect and early discrimination, thereby promoting the 
emergence of pre-mating barriers among flycatcher popula- 
tions.?^/ Although perceptual predispositions expressed early in 
life can be significantly altered by experience, "^^? weak re- 
sponses to songs outside the typical range of the local dialect, 
irrespective of species, imply that dialect-specific responses 
are likely to promote incipient pre-mating reproductive barriers 
between the collared and pied flycatcher prior to secondary 
contact. 

Previous work on collared flycatchers identified the earliest 
metabolic responses to songs around day 4, while species 
discrimination did not develop until around day 9.°° This is 
consistent with what is known about auditory development in 
the pied flycatcher. Electrophysiology experiments on nestling 
birds demonstrated responses in the auditory cortex to tone 
frequencies typical of pied flycatcher vocalizations between 
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Figure 3. Playback categories and nestling responses to local- 
similar song playbacks 

(A) Euclidean distance between playbacks and local dialect centroids. Light 
gray, local songs; dark gray, foreign songs; blue, collared flycatcher songs. 
Playback categories denoted by dashed lines. Local-similar category contains 
9596 of local dialect songs, while collared-similar category contains 9096 of 
collared flycatcher songs. 

(B) Responses of 11-day-old nestlings to the different playback categories. 
(C) Responses of 12-day-old nestlings to the different playback categories. 
Clutches that produced no begging calls during any playback were removed 
for plotting purposes. Responses are averaged across nestlings within 
clutches. Boxplots give median, quartiles, and quartiles + 1.5 x IQR. p values 
are from post hoc tests (see Figure S1 and Table S3). 


days 4 and 5,°°°” whereas specific responses to natural pied 
flycatcher vocalizations did not arise until day 11.°° Thus, the 
nestling auditory system first attains the ability to perceive 
songs by day 4, whereas discrimination among behaviorally 
relevant sounds develops between days 9 and 11. Our findings 
that song discrimination increases with nestling age are consis- 
tent with these earlier studies. The apparent earlier expression 
of species level (by day 11) compared with dialect-level 
discrimination (by day 12) implies that, like other songbirds, 
the pied flycatcher auditory system achieves increasingly spe- 
cific responsiveness with age.?? Next, we explored whether 
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increasing specificity is likely to be a result of early exposure 
to conspecific sounds. 


No evidence for an effect of song experience on song 
responses 

The responses of nestling songbirds to the songs of heterospe- 
cifics have been shown to be either independent of exposure to 
the other species’ songs ^^ or independent of the early environment 
in general.'* The degree to which dialect-specific responses 
depend on previous experience is less clear. Nestling golden- 
crowned sparrows express stronger responses to their local dia- 
lect,”° but experiments on the closely related white-crowned spar- 
row suggest that dialect-specific responses depend on previous 
experience.?? Thus, nestling song responses to playbacks similar 
to the local dialect might arise due to learning from early song 
exposure from surrounding males, including the social father. 

We used three tests to evaluate whether nestling song re- 
sponses are affected by early exposure to conspecific songs. 
First, in a single population, Tovetorp, we estimated early song 
exposure using automated recordings at 28 nest boxes during 
the nestling period (https://doi.org/10.17632/v49k6v4fw3.1). 
We found weak evidence for an effect of song exposure on nest- 
ling song responses (likelinood-ratio test: x2 = 3.61, p = 0.165). 
Contrary to the prediction that song exposure to local songs 
should increase discrimination, post hoc tests suggested that 
nestling responses to local-similar songs decrease with 
increasing song exposure (Table S4). Second, at 24 of these 
nest boxes, we measured the acoustic features of songs present 
in these recordings, presumed to be from the nestlings' social fa- 
ther (https://doi.org/10.17632/v49k6v4fw3.1). Contrary to the 
prediction of early learning promoting song responses, we found 
evidence that nestlings respond weaker to playbacks similar to 
the songs of their social father (estimated coefficient on dissim- 
ilarity to play back: 5.13 + 2.82 SE, likelihood-ratio test: x; =3.19, 
p = 0.074; Table S4). Finally, we tested whether the similarity of 
the social father’s song to the local dialect influences responses 
to local-similar songs. We found little evidence for an effect of the 
similarity of the father's song to the local dialect (likelihood-ratio 
test: X2 = 0.52, p = 0.771). Again, contrary to the prediction of 
early learning promoting song responses, post hoc tests sug- 
gested trends for nestling responses to local-similar songs to 
decrease with the similarity of the father's song to the local dia- 
lect (Table S4). 

We found no evidence that nestling responses to local dialect 
songs are increased by early song experience. In contrast, our 
results suggest that nestlings respond less to playbacks similar 
to their social fathers' songs. The context and function of male 
singing during the nestling stage are unclear, but one possibility 
is that males sometimes produce song before feeding visits, ^' 
leading nestlings to naturally produce begging calls when hear- 
ing their social father's song. Males that sing at abnormally high 
rates during the nestling phase —for example, due to seeking out 
extra-pair copulations? —may decouple song production from 
feeding and thereby drive habituation of the begging response 
in highly exposed nestlings. This does not imply that nestlings 
do not learn songs from their father.*' Rather, it suggests that 
early responses to dialect variation are unlikely to be explained 
through early learning. Without experimental manipulation of 
early experience, we remain cautious in our conclusions. 
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However, no matter the mechanism driving it, nestlings express 
early discrimination in favor of their own dialect shortly before 
leaving the nest and well before song learning begins in pied 
flycatchers.? 


Conclusion 

We show that pied flycatcher nestlings respond most 
strongly to songs with acoustic qualities typical of their own 
local dialect, an ability that emerges just before fledging. We 
demonstrate that co-occurrence with the collared flycatcher 
has little effect on early song responses in pied flycatchers, 
suggesting that the risks of cross-species mating have, at 
most, a minor and secondary influence on the evolution of 
early song discrimination. Instead, we show that nestling birds 
express reduced responses to songs that are dissimilar 
from their local dialect, irrespective of from which species 
they derive. This implies that species discrimination is a 
pattern resulting from acoustic dissimilarity of songs, 
rather than the costs of species interactions. Although the 
precise developmental mechanisms leading to associations 
between songs and early discrimination within populations 
remain to be clarified, the negligible effects of early experi- 
ence on nestling song responses suggest innate develop- 
ment. Our results therefore support the long-held hypothesis 
that cultural evolution can drive the emergence of song 
learning predispositions, which hastens the divergence in 
mate recognition systems between closely related species.^? 
The tremendous species richness of vocal learning birds ^? 
may thus be partially explained as a side-effect of cultural 
evolution. 
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EXPERIMENTAL MODEL AND SUBJECT DETAILS 


Pied flycatchers are small, insectivorous, migratory songbirds.? They winter in central and western Africa, while their breeding range 
in Europe extends from Spain and the United Kingdom in the west to northern Scandinavia and to central Russia in the east. Song 
playback experiments were conducted between 2018-2021 at six distinct breeding populations along a latitudinal gradient in the 
western portion of the pied flycatcher's range: Valsaín (Spain V, 40°52’N, 4°01’W, 2019 and 2020), La Hiruela (Spain LH, 41°4’N, 
3°27’W, 2018 and 2020), East Dartmoor (England, 50°36’N, 3°43’W, 2018), De Hoge Veluwe National Park (The Netherlands, 
52°04’N, 5°49’E, 2018), Tovetorp (Sweden, 58°57’N, 17°089’E, 2020 and 2021), and Turku (Finland, 60°27’N, 22°15’E, 2018). In addi- 
tion, song playback experiments, part of a previous study, ^ were conducted between 2013-2016 on Öland, Sweden (57°10°N, 
17*00'E). Geographic distances among playbacks and populations were calculated using the R package geodist 0.0.4.^/ The closest 
sympatric site to Tovetorp and Finland was Oland, while the closest sympatric site to all other populations was estimated in north- 
eastern France (48°56’N, 5°1’E). 
Animal procedures were approved by Jordbruksverket (Linkópings djurfórsóksetiska nàmnd, Dnr 01110-2020). 


METHOD DETAILS 


Song recordings and playbacks 

Together, song recordings from 168 different individuals were analyzed. Each recording contained songs from a single individual. 
Most individuals were recorded from the populations where subsequent playback experiments were conducted. Four individuals 
from the Netherlands, one individual from England, 7 collared flycatchers, and all of the songs from Finland were obtained from 
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an online repository of avian vocalizations (https://xeno-canto.org/). For all other populations, recordings were made by the authors 
with a condenser microphone (Sennheiser K6 module and ME66 capsule, Sennheiser Electronic GmbH & Co. KG, Wedemark, Ger- 
many) connected via XLR cable to a digital audio recorder (Tascam DR-40, TEAC Corporation, Tokyo, Japan or Fostex FR-2, Foster 
Electric Company, Ltd, Tokyo, Japan). For each pied flycatcher population, we obtained 11-30 recordings (median: 15), whereas we 
obtained 50 recordings of collared flycatcher individuals from multiple populations. 


Song analysis 

We analyzed between 2-34 song phrases for each of the 168 individuals for which we had song recordings (median = 9). Song re- 
sponses are likely to be determined by these acoustic features, rather than species identity or population identity per se. In order to 
extract a measure of song variation independent of these categories, we analyzed and compared songs from both species and from 
each pied flycatcher population. Songs were measured by a single author (D.W.) using Luscinia software. ^^ Each song recording was 
imported and song phrases were measured separately. Spectrograms were visualized using a Gaussian windowing function with the 
following settings: 12 KHz maximum frequency, 5 ms frame length, 220 spectrograph points, 80% spectrograph overlap, 50 dB dy- 
namic range, 3096 dereverberation, and 50 ms of dereverberation range. For each song phrase, we measured individual elements as 
continuous sound traces. Elements were grouped into syllables if their sound traces were overlapping in time or the gap between 
them was less than 30 ms. 

We compared songs using Luscinia's built-in dynamic time-warping (DTW) algorithm, which finds the optimal alignment of each pair 
of syllables based on multiple acoustic features and uses this alignment to calculate dissimilarity.^^ Songs are subsequently compared 
by aligning their syllables sequentially. We performed DTW using the following acoustic feature weightings: time (5), mean frequency 
(1), mean frequency change (1), and normalized mean frequency (1). All other measures were weighted by their standard deviation. We 
used a compression factor of 0.001, a minimum element length of 25 samples, time SD weighting of 1, a maximum warp of 2596, and 
compared syllables by stitching elements. All other options were left at the default settings. These settings were chosen based on the 
recommended values and because they reliably grouped song phrases from the same individual that were subjectively determined to 
be highly similar. Luscinia generates dissimilarity matrices for each acoustic measure and then performs non-metric multi-dimen- 
sional scaling to reduce the number of dimensions while preserving overall dissimilarity. Subsequently, Luscinia performs a principal 
components analysis to generate 10 independent axes of variation. The eigenvalues of each principal component relate to the pro- 
portion of total variance in dissimilarity explained. Ten principal components explained 75% of the total variance in dissimilarity among 
song phrases. A Kruskal stress test had a stress value of 0.025 indicating a good representation of dissimilarity among songs. 

Due to the potential inaccuracy in acoustic measurements of recordings in mp3 format (i.e., those taken from https://xeno-canto. 
org/), we repeated all analyses excluding these recordings and found highly similar results. Thus, all recordings were included in the 
results presented. 

The 10 principal component scores were used to compare playbacks (taken as the mean of the component songs within a play- 
back), individuals (taken as the mean centroid of the songs within an individual), population (taken as the mean centroid of the indi- 
viduals within the population), and species (taken as the mean centroid of the individuals within the population). Songs varied 
between species in PC1 (linear mixed models with individual nested within population as random effects; N = 2088 song phrases, 
N = 168 individuals, x; = 11.58, P « 0.001) and there was a weak trend for species-level differences in PC2 (x? = 1.98, 
P = 0.159), but we did not detect evidence for differences in the other principal components (likelihood-ratio test: x? « 1.66, 
P > 0.198). PC1 scores alone correctly classified 100% of the collared flycatcher singers (n = 41 individuals, mean PC1 
score + SD, 0.11 + 0.02) and 97.6% of the pied flycatcher singers (n = 127, -0.02 + 0.04) to the species level. Given this independent 
support for species-level song classification, we considered playback species as a categorical variable in subsequent analyses. 

We found that songs varied across populations in 9 out of 10 principal components (linear mixed models with individual as random 
effect for PC1-PCS, likelihood-ratio test: X > 14.20, P < 0.015). We calculated the Euclidean distance of each playback in multidimen- 
sional principal component space to the mean centroid of each population and used this as an estimate of a playback’s similarity to 
songs from a given population. Unsurprisingly, pied flycatcher playbacks are closer to the centroid of their source population (n = 80 
playbacks; mean Euclidean distance + SE: 0.06 + 0.00) than to the centroids of other populations (n = 143; 0.10 x 0.00): 95% of all 
playbacks are within 0.10 of their source population’s centroid, while around 50% of foreign playbacks are within 0.10 of other pop- 
ulations' centroids. Likewise, collared flycatcher playbacks have greater Euclidean distances to population centroids compared to 
either local or foreign playbacks (n = 31; 0.17 + 0.00; Figure S1B): 90% of collared playbacks, but only 3% of pied flycatcher playbacks, 
are more than 0.15 from population centroids. We considered nestling responses to playback similarity to songs from their own pop- 
ulation using Euclidean distance as a factor (‘“‘local-similar”: < 0.10, *local-dissimilar": > 0.10 and < 0.15, "collared-similar": > 0.15). 


Playback recordings 

One-minute-long playback files were obtained directly from the song recordings used to compare species and populations. Individ- 
ual recordings that were contained more than three minutes of song were divided up into distinct playbacks. In total, we produced 
245 playbacks from 168 individuals. A single author (D.W.) removed background sounds using RAVEN v1.5“° and constructed play- 
back files containing between 8-12 song phrases, each separated by four seconds of background noise. Background noise was 
obtained either from the recording itself or from recordings from the same site, in which case the amplitude was adjusted to match 
the background of the song recording. At each site, we played back between 26-31 unique collared flycatcher playbacks (median 
30.5), 27-75 unique foreign playbacks (38.5), and 12-36 unique local playbacks (20.5). 
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Playback experiments 

We broadcast pied and collared flycatcher song playbacks to whole clutches of pied flycatcher nestlings (N = 347 clutches, for a total 
of 1722 nestlings). Each nestling was marked with a unique symbol on the top of its head using correction fluid to allow individual 
discrimination during video analysis and placed into an experimental wooden nest box. The nest box had an artificial nest con- 
structed of moss and/or leaves on the bottom and was fitted with a video camera (ELP-USB-130W01MT-MGAO[BW], Alipu Technol- 
ogy Co., Ltd, Guangdong, China) and a microphone (Go Mic, Samson Technologies, Hicksville, NY, USA) positioned to record the 
behavior and vocalizations of all nestlings. The camera and microphone were connected to a laptop computer for recording. A 
speaker was positioned roughly one meter away and facing the nest box for song playbacks (Eco Extreme, Grace Digital, Poway, 
CA, USA). Around half of the experiments conducted at Tovetorp used a different video recording device, as described in Wheatcroft 
and Qvarnström. '^ 

Nestlings acclimatized to the nest box for at least 5 minutes before the experiment started. An experiment consistent of three one- 
minute periods of silence each followed by a one-minute song playback. The order of the song playback treatments (collared, foreign 
and local pied flycatcher) was alternated between experiments. The individual song playbacks making up an experiment were alter- 
nated to limit repeated pairing of the same playbacks. Throughout each experiment, audio and video of the nestling behaviors were 
recorded. Song playbacks were broadcast using a smartphone connected to the speaker. Playbacks were broadcast at 75 dB, 
measured at 3 m distance (Digital Sound Level Meter, RadioShack). After the third playback, the nestlings were weighed and returned 
to their home nest box. Most (n = 327 out of 347) clutches were tested once. Broods were tested on a second occasion, at least 12 
hours after the first, only when no nestlings produced any behavioral response on the first trial (Tovetorp: n = 20 clutches repeated). 
Between 29 (Finland) and 99 (Spain V) clutches were tested at each population (median: 38, supplemental information). In addition, 
37 clutches were tested on Óland. 

We recorded individual mass measurements around the time of playback for most individual nestlings (n = 1718 out of 1722 nes- 
tlings, excluding Oland). In order to include behavioral responses for those nestlings without individual mass measurements, we 
imputed their mass by taking the mean mass of their nestmates. For nestlings with individual mass measurements, a large proportion 
of the variance in individual mass was explained by the mean mass of nestlings in the clutch (adjusted /? in a linear model = 0.73, 
estimated coefficient on mean mass = 1.00 + 0.02 SE). For nestlings from Oland (n = 186 nestlings), where we could not associate 
individual mass with their nestling behavioral responses, we used the mean nestling mass of each clutch. Two experiments from To- 
vetorp, one from Spain LH, and seven from Óland were excluded from subsequent analyses due to missing nestling mass data. 


Behavior analysis 

A single observer (L.B.) watched and noted the behavioral responses of all nestlings in each video recording both during periods of 
silence and song playbacks. Individual nestlings were tracked by the painted symbols marked on their heads. These symbols could 
later be associated with ring numbers and masses recorded on the day of playback. The observer was partially blind to the exper- 
imental treatment based on initial unfamiliarity with distinguishing the songs of either species and inability to distinguish the songs by 
population, but could identify when songs were being played back as well as the population from which each video originated. Based 
on previous work in this system, '^?* the observer recorded the following behaviors as counts during each silence and song playback 
period: 1) begging calls (defined as an audible low-intensity begging call, Video S1), 2) gapes (defined as opening and closing the bill 
completely), 3) jumping begging calls (defined as a high-intensity begging call, accompanied by jumping forward or upward), 4) po- 
sition shifts (defined as the nestling moving continuously from one orientation or position to another, such that there was little overlap 
in its body after shifting), and 5) looking up (defined as the nestling shifting the orientation of its head and bill from perpendicular to 
facing the nest box opening). For begging calls and jumping begging calls, the observer referred to spectrograms of extracted audio 
to confirm counts when multiple nestlings were begging at the same time. 

Producing begging calls and looking up were the most common behaviors (n = 540 and 375 out of 1908 nestlings produced at least 
one begging call or look-up, respectively, during one of the playback treatments, including the Oland population). During a given play- 
back treatment, these two behavioral measures were positively correlated within nestlings (N = 5516 playback treatments, Spear- 
man's rho = 0.47, P « 0.001). We used begging calls as a response variable in subsequent analyses to facilitate comparison with 
previous studies and because it was determined to be less subjective than looking up. ^^? 

Videos of pied flycatcher nestlings from Oland (n = 37 clutches, for a total of 186 nestlings), collected using equivalent methods to 
the current study, were re-analyzed by L.B. These data were originally published in Wheatcroft and Qvarnstrém.'* We conducted 
tests of repeatability of begging call counts made by L.B. and Wheatcroft and Qvarnstrém. '“ Since nestlings in the Öland population 
were not individually marked, we summarized counts within each trial by taking the mean of nestling counts within each trial. We as- 
sessed repeatability of begging call counts across observers using rptR v0.9.2.^? Begging call counts were highly repeatable across 
observers for both begging calls prior to (1000 bootstrap iterations, R = 0.88, 95% Cl 0.83 — 0.92, P < 0.001) and during the playback 
period (1000 bootstrap iterations, R = 0.89, 95% CI 0.83 - 0.92, P < 0.001) 


Estimating song exposure 

Audio recordings were conducted at Tovetorp in 2020 and 2021 for around three hours in the morning (approximately two hours 
around civil dawn and one hour between 7:00 - 8:00AM) and we obtained multiple recordings for each nest box (12.5 median; range: 
4.4 - 22.3 hours of recordings per box). A single observer (D.W.) annotated spectrograms generated from 10 two-hour-long record- 
ings in RAVEN,*° forming selection boxes enclosing each song phrase in time and frequency. This resulted in 1,108 song phrase 
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selections from 10 individual males. Standardized images of each selection (1000 x 500 pixels, 10 s x 15 kHz) were generated using 
custom scripts in Python and used to train (8096 of images), validate (1096), and test (1096) a cloud-based object detection model 
(https://cloud.google.com/vertex-ai/). The model had an average precision of 0.903. Evaluation of the false negatives suggested 
that the model tended to miss song phrases that had just begun or were soon to end within the image. Overlapping windows 
were used for model predictions to increase the detection rate. 

The model was downloaded and run within a docker container (https://www.docker.com/). Standardized spectrograms from 
each audio recording were generated (500 x 500 pixels, 5 s x 15 kHz) with a 2.5 s window overlap. We obtained model pre- 
dictions for each image, using custom Python scripts to account for window overlap. Predictions were filtered with a confidence 
score of 0.95 and all predictions were subsequently manually validated through analysis of audio and/or spectrogram. Over all 
boxes, we recorded an average of around one song every 10 minutes, with substantial variation across nests (0.08 + 0.19 
mean songs per minute x SD). We estimated song exposure as the total number of song phrases divided by the total recording 
time. Counts of song phrases were divided by the total recording time to obtain estimated song exposure. To obtain acoustic 
features of the songs to which each clutch was exposed, we extracted song phrases from audio recorded at 24 nest boxes 
(median: 9.5 song phrases per nest box, range: 1 - 16) and analyzed the songs in Luscinia using identical methods as above. 
All songs were subject to a second dynamic time warping procedure with identical methods as above, and Euclidian distances 
between the songs to which nestlings were exposed and playbacks were calculated using all 10 principal components (see 
Table S4). 


QUANTIFICATION AND STATISTICAL ANALYSIS 


All statistical tests were performed in R v3.6.1.^? To evaluate differences in song principal component scores between species and 
among populations, we used linear mixed models, implemented using /mer in Ime4 v1.1-21,^? using individual singer as a random 
intercept and species or population as a fixed effect. Significance of interactions and main effects was evaluated using likelihood 
ratio tests. 

We evaluated the correlation between genetic, geographic, and song divergence among populations using the function 
MMRR.°° 

To evaluate variation in begging calls in response to songs, we first tested for overdispersion of begging calls using the function 
dispersiontest in package AER v. 1.2-7°' on fixed effect poisson generalized linear mixed models with begging calls during the play- 
back as the dependent variable. We found evidence for significant overdispersion. Similarly, visualization of the data revealed that a 
large number of trials consisted of zero responses (i.e., no nestling in a given clutch produced a begging call during a given playback 
treatment). We then applied zero-inflated quasi-poisson regression using the R package glmmTMB v1.0.2.1.°7 In all models, we 
included the order of the given playback within a trial, the brood size, estimated nestling mass within a given nest (g), the number 
of begging calls produced by each individual nestling during the immediately preceding silent period, and the age of the nestlings. 
In addition, we included a random effect of the playback recording used to stimulate begging calls: 


(a) begSpiayback ~ intercept + order + brood. size + mass + begSsilence + age + (1|playback) 


In addition, all models included a zero-inflation term in each model that consisted of an intercept and a random effect of the brood 
identity: 


(b) begSpiayback ~ intercept + (1|brood_id) 


In all subsequent models, we evaluated interactions between age and relevant song measures to estimate the ontogeny of discrim- 
ination (see Table S1). To test song discrimination in allopatry, we included the playback species as a fixed effect in the above model. 
To compare song discrimination in sympatry and allopatry, we included an interaction term between the playback species and either 
sympatry, as a categorical variable, or distance to sympatry, as a continuous variable (see Table S2). To test the influence of playback 
similarity to the local dialect, we included the similarity category (*local-similar", *local-dissimilar", *collared-similar") as a fixed ef- 
fect (see Table S3). Finally, to test the effect of song exposure on song responses at Tovetorp, we included interactions between song 
exposure, estimated as the total number of song phrases divided by the total recording time, and the Euclidian distance between the 
playback and the social father's song as fixed effects (see Table S4). Due to convergence errors, we did not include random effects 
for models including only Tovetorp. 

For all models, we tested the significance of fixed effects and interactions using log-likelihood ratio tests and, to evaluate interac- 
tions, we performed post-hoc tests using the functions emmeans and emtrends in the emmeans package v1.4.3.01.°° For interac- 
tions between two categorical variables, we evaluated the difference between levels of one variable across levels of the other variable 
using the contrast function in emmeans. For interactions between two continuous variables, we evaluated the effect of one variable at 
set levels of the other variable (mean, mean x SD). 

The majority of experiments (8496 of N 2 347) were conducted at nests with nestlings age 11 or 12, but around twice as many ex- 
periments were conducted on age 12 (n = 195) compared to age 11 nestlings (n = 98). In addition, nests with age 11 nestlings were 
moderately less likely to produce begging calls (4896 of nests) compared to nests with age 12 nestlings (5996 of nests). As a result, 
the observed increase in song discrimination with age (Figures 3B and 3C; Table S3) could be partially explained by sample size 


e4 Current Biology 32, 1-6.e1—-e5, December 5, 2022 


Please cite this article in press as: Wheatcroft et al., Species-specific song responses emerge as a by-product of tuning to the local dialect, Current 


Biology (2022), https://doi.org/10.1016/j.cub.2022.09.063 


Current Biology © Cell?ress 


OPEN ACCESS 


and begging rate differences across age categories. To evaluate this possibility, we bootstrapped our findings by generating 100 
random sub-samples of experiments at nests with age 12 nestlings. Each sub-sample consisted of 80 nests ([begging rate at age 
11 nests]/[begging rate at age 12 nests] x [sample size of age 11 nests]), randomly selected with replacement. We applied identical 
zero-inflated quasi-poisson regression models to all 100 datasets and evaluated the difference in song discrimination between age 
11 and age 12 as for the entire dataset. The effect of age on song discrimination was evaluated as the proportion of model runs 
with a significant difference in discrimination between ages 11 and 12. 
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Figure S1. Overview of main methods, Related to Figures 1, 2, and 3, and STAR Methods. 
(A) Acoustic analysis on song phrases used to generate principal component scores for 
playbacks, individuals, populations, and species. (B) Playback experiments. 


Estimated coefficient + SE Z P 


Conditional model 


Intercept 1.28 + 0.86 1.50 0.134 
Beg.rest 0.04 + 0.00 16.56 « 0.001 
Estimated.mass -0.16 + 0.03 -6.16 « 0.001 
Playback.order -0.14 + 0.05 -2.89 0.004 
Clutch.size -0.06 + 0.03 -1.77 0.076 
Age 0.25 + 0.07 3.72 « 0.001 
Playback.speciesPied 0.33 + 0.10 -3.46 « 0.001 
Post-hoc tests t ratio 

Agell:Pied 0.29 € 0.15 1.96 0.051 
Agel2:Pied 0.33 + 0.10 3.48 « 0.001 
Agel3:Pied 0.37 + 0.15 2.43 0.015 


Table S1. Nestling begging calls explained by playback, Related to Figure 2 and STAR 
Methods. Results of a zero-inflated quasi-poisson mixed effect model with nestling begging call 
as dependent variable and playback species as fixed effect. Beg.rest refers to the begging calls 
produced by the nestling during the minute preceding the song playback. Post-hoc tests 
(estimated marginal means) estimate the effect of age on responses to playback species. 


Estimated 
coefficient + SE t ratio F 
Sympatry distance (km) 


Dist86:Pied 0.44 + 0.15 2.91 0.004 
Dist435:Pied 0.38 + 0.10 3.38 < 0.001 
Dist784:Pied 0.32 + 0.10 3.17 0.002 
Sympatry 

Allopatry 0.32 + 0.10 3.20 0.001 
Sympatry 0.60 + 0.40 1.49 0.136 


Table S2: Effect of sympatry on nestling responses to playback species, Related to Figure 2. 
Post-hoc tests assessing effect of playback species on nestling responses based on population 
distance to sympatry as well as in sympatry as a factor using estimated marginal means. 


Estimated coefficient SE Zz P 


Conditional model 


Intercept 1.77+ 0.87 2.03 0.042 
Beg.rest 0.04 + 0.00 16.44 < 0.001 
Estimated.mass -0.16 + 0.03 -6.22 < 0.001 
Playback.order -0.15 + 0.05 -3.20 0.001 
Clutch.size -0.06 + 0.03 -1.80 0.073 
Age 0.25 + 0.07 3.73 « 0.001 
Category -2.05 + 0.96 -2.14 0.032 
Post-hoc tests tratio P 
LD:CS 0.02 + 0.11 0.18 0.983 
LS:CS 0.24 € 0.11 2.14 0.081 
LS:LD 0.22 + 0.09 2.42 0.041 
Post-hoc tests tratio P 
Agell 

LD:CS 0.40 + 0.18 2.24 0.065 
LS:CS 0.25 + 0.17 1.41 0.337 
LS:LD -0.15 + 0.16 -0.96 0.603 
Agel2 

LD:CS 0.01 +0.11 0.04 0.999 
LS:CS 0.23 +0.11 2.09 0.092 
LS:LD 0.23 + 0.10 2.48 0.035 
Agel3 

LD:CS -0.39 + 0.19 1.99 0.114 
LS:CS 0.22 + 0.18 1.20 0.456 
LS:LD 0.60 + 0.16 3.73 0.001 


Table $3: Nestling begging calls explained by playback similarity to local dialect, Related to 
Figure 3 and STAR Methods. Results of zero-inflated quasi-poisson mixed effect models with 
nestling begging call variation as dependent variable and local dialect similarity as fixed effect. 
Post-hoc tests estimating responses to playbacks with varying local dialect similarity 
(“Category”) and the effect of nestling age on these responses using estimated marginal means. 
“LD”: local-dissimilar, “LS”: local-similar, “CS”: collared-similar. 


Model 1: effect of song exposure on song responses 


Post-hoc tests Estimated contrast + SE t ratio T 

LS -9.60 + 7.36 -1.31 0.193 
LD -1.89 + 1.67 -1.13 0.258 
CS -0.54 € 1.12 -0.57 0.571 


Model 2: effect of similarity between social father's song and playback 


Estimated coefficient + SE Z value PD 
Intercept 3.69 + 3.90 0.95 0.344 
Beg.rest 0.07 + 0.01 10.42 « 0.001 
Estimated.mass -0.02 + 0.16 -0.15 0.878 
Playback.order -0.30 + 0.19 -1.61 0.108 
Clutch.size -0.12 + 0.12 -1.01 0.314 
Age -0.22 + 0.24 -0.90 0.368 
Father.ED.to.Playback 5.13 + 2.82 1.82 0.069 


Model 3: effect of similarity between social father's song and local dialect on song 


responses 
Post-hoc tests Estimated contrast + SE t ratio P: 
LS 6.45 + 5.60 1.15 0.251 
LD 1.29 € 7.32 0.176 0.860 
CS 1.54 + 5.61 0.275 0.784 


Table S4: Nestling begging call responses explained by the social father's song, Related to 
STAR Methods. Results of zero-inflated quasi-poisson models with nestling begging call 
variation as dependent variable. Model 1: post-hoc tests estimating the effect of varying levels of 
song exposure on nestling responses to the playback similarity to the local dialect. Model 2: 
model coefficients for a model estimating the effect of the Euclidean distance of the social 
father's song to the playback. Model 3: post-hoc tests estimating the effect of the similarity of 
the father's song to the local dialect. 


